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Abstract 

We present a recent analysis of e'/e in the 1/N C expansion. We show 
that the 1/N C corrections to the matrix element of Qa are large and 
positive, indicating a AI — 1/2 enhancement similar to the one of Qi 
and Q2 which dominate the CP conserving amplitude. This enhances 
the CP ratio and can bring the standard model prediction close to the 
measured value for central values of the parameters. 
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We present a recent analysis of e'/e in the 1/N C expansion. We show that the 
1/N C corrections to the matrix element of Qg are large and positive, indicating 
a AI = 1/2 enhancement similar to the one of Qi and Q2 which dominate the 
CP conserving amplitude. This enhances the CP ratio and can bring the standard 
model prediction close to the measured value for central values of the parameters. 



1 Introduction 



Direct CP violation in JC — > tttt decays was recently jpbserved by the KTeV and 
NA48 collaborationsBEl The present world averageB for the parameter e'/e is 
Re e'/e — (21.2 ± 4.6) • 10~ 4 . In the standard model CP violation originates 
in the CKM phase, and direct CP violation is governed by loop diagrams of 
the penguin type. The main source of uncertainty in the calculation of e'/e is 
the QCD non-perturbative contribution related to the hadronic nature of the 
K — > tttt decay. Using the AS = 1 effective hamiltonian, 



<i/AS=l 

H eff 



71 u 



t=i 



< m c ) , 



(1) 



the non-perturbative contribution, contained in the hadronic matrix elements 
of the four-quark operators Qi, can be separated from the perturbative Wil- 
son coefficients Cj(/x) = Zi(fi) + rj/j(/z) (with r = — A t /A„ and X q = V* s V qd ). 
Introducing (Qi)/ = ((7T7r)j|Qi|.K'), the CP ratio can be written as 



Gf u> lmXt 
~Y \e\ReA 



(2) 



u) =RcAo/Rcyl.2 = 22.2 is the ratio of the CP conserving K — > tttt isospin arru 
plitudes; f2^+^' encodes the effect of the isospin breaking in the quark massesB 
e'/e is dominated, by (Qq)q and {Qs)2 which cannot be fixed from the CP 
conserving dataQ^El Beside the theoretical uncertainties coming from the calcu- 
lation of the {Qi)i and of £1^+^ , the analysis of the CP ratio suffers from the 
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uncertainties on the values of various input parameters, in particular of the 
CKM phase in ImAt, of A QCD = A^i, and of the strange quark mass. 

To calculate the hadronic matrix elements we start from the effective chi- 
ral lagrangian for pseudoscalar mesons which involves an expansion in mo- 
menta where terms up to 0(p 4 ) are includedH The method we use is the l/N c 
expansionQu In this approach, we expand the matrix elements in powers of 
the momenta and of XlMc- For the l/N c corrections we calculated chiral loops 
as described in refs. BlJ. Especially important to this analysis are the non- 
factorizable corrections, which are UV divergent and must be matched to the 
short-distance part. They are regularized by a finite cutoff A c which is iden- 
tified with the short-distance renormalization scale. The definition of the mo- 
menta in the loop diagrams, which are not momentum translation invariant, is 
discussed in detail in ref . u. Other recent work on matrix elements in the 1 /N c 
approach can be found in refs.EiEj. 

For the Wilson coefficients we use the leading logarithmic and the next-to- 
leading logarithmic valuesu The absence of any reference to the renormaliza- 
tion scheme in the low-energy calculation, at this stage, prevents a complete 
matching at the next-to-leading ordert^l Nevertheless, a comparison of the nu- 
merical results obtained from the LO and NLO coefficients is useful as regards 
estimating the uncertainties and testing the validity of perturbation theory. 

2 Analysis of e' / e 

Analytical formulas for all matrix elements, at next-to-leading order in—the 
twofold expansion in powers of momenta and of 1/N C , are given in refs.Hil3. 
In the pseudoscalar approximation, the matching has to be done below 1 GeV. 
Varying A c between 600 and 900 MeV, the bag factors B^ 2 ^ and B^ 2 ^ take 
the values 8.2— 14.2 and 2.9 — 4.6; quadratic terms in (Qi)q and (Q2)o produce 
a large enhancement which brings the AI — 1/2 amplitude in agreement with 
the dataEHl Corrections beyond the chiral limit were found to be small. 

For {Qe)o and (Qg)2 the leading non-factorizable loop corrections, which 
are of O(p /N c ), are only logarithmically divergently Including terms of O(p ), 
0{p 2 ), and O(p /N c ), B { 6 1/2) and Bf /2) take the values 1.10 - 0.72 and 
0.64— 0.42. As a result the experimental range for e'/e can be accommodated 
in the standard model only if there is a conspiracy of the input parameters. | 
However, since the leading 0(p ) contribution vanishes for Qq, corrections from 
higher xu-der terms beyond the 0(p 2 ) and O(p /N c ) are expected to be large. 
In ref. t3 we investigated the 0(p 2 /N c ) contribution, i.e., the 1/N C correction 
at the next order in the chiral expansion, because it brings about, for the 
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Table 1: Numerical values for e'/e (in units of 10 4 ) as explained in the text. 



LO 


14.8 


< 


e'/e 


< 


19.4 


6.1 


< 


e'/e 


< 


48.5 


NDR 


12.5 


< 


e'/e 


< 


18.3 


5.2 


< 


e'/e 


< 


49.8 


HV 


7.0 


< 


e'/e 


< 


14.9 


2.2 


< 


e'/e 


< 


38.5 





















first time, quadratic corrections on the cutoff. From counting arguments and 
more generally from the fact that the chiral limit is assumed to be reliable, the 
quadratic terms (which are not chirally suppressed) are expected to be domi- 
nant. It is still desirable to check that explicitly by calculating the corrections 
beyond the chiral limit, from logarithms and finite terms, as done for Q\ and 
Qi. Numerically, we observe a large positive correction -from the quadratic 
term in (Qe}o- This point was already emphasized in ref. t3. The slope of the 
correction is qualitatively consistent and welcome since it compensates for the 
logarithmic decrease at O(p /N c ). Varying A c between 600 and 900 MeV, the 
B^ 2 ^ factor takes the values 1.50— 1.62. Qq is a AI = 1/2 operator, and the 
enhancement of {Qe)o indicates that at the level of the 1/N C corrections the 
dynamics of the AI =1/2 rule applies to Qq as to Q\ and Qi. 

Using the quoted values for B^ 2 ^ 1 together with theiiill leading plus next- 
to-leading order B factors for the remaining operatorsEJ we calculated e'/e. 
The results for the three sets of Wilson coefficients LO, NDR, and HV and for 
A c between 600 and 900 MeV are given in Tab. |[ The numbers are close to the 
measured value for central values of the parameters (first column). They are 
obtained by assuming zero phases from final state interactions. This approxi- 
mation is very close to the results we, would get if we used the small imaginary 
part obtained at the one-loop levelE-3 

Performing a scanning of the parameters [125 MeV < m s (lGeV) < 175 
MeV, 0.15 < Q. v+r) , < 0.35, 1.04 • 10" 4 < ImA t < 1.63 • 10" 4 , and 245 MeV < 
Aqcd _ 405 MeV] we obtain the numbers in the second column of Tab. [l]. They 
can be compared with the results of refs.!'0llllll'EJ The values of Bg 1 ^ 2 ^ can 
also be compared with ref. and those of B { 8 3/2) with refs. 00. The large 
ranges reported in the table can be traced back to the large ranges of the input 
parameters. This can be seen by comparing them with the relatively narrow 
ranges obtained for central values of the parameters. The parameters, to a 
large extent, act multiplicatively, and the large range for e'/e is due to the 
fact that the central value(s) for the ratio are enhanced roughly by a factor 
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Table 2: Same as in Tab. |l[ but now with the phenomenological values for the phases as 
explained in the text. 



LO 


19.5 


< 


e'/e 


< 


24.7 


8.0 


< 


s 




< 


62.1 


NDR 


16.1 


< 


e'/e 


< 


23.4 


6.8 


< 


e 


h 


< 


63.9 


HV 


9.3 


< 


e'/e 


< 


19.3 


2.8 


< 


e 




< 


49.8 























of two compared to the results obtained with B factors for Q$ and Qg close 
to the VSA. More accurate information on the parameters, from theory and 
experiment, will restrict the values for e'/e. 

To estimate the uncertainties due to higher order final state interactions 
we also calculated e'/e using the real— part of the matrix elements and the 
phenomenological values of the phasesE3, <5 = (34.2 ± 2.2)° and 82 — (— 6.9± 
0.2)°, i.e., we replaced |£,-W<Qi>/l in Eq. @ by £. yiRe(Qi) z / cos jj. The 
corresponding results are given in Tab. §. They are enhanced by - 25 % 
compared to the numbers in Tab. [l]. We would like to emphasize that this 
~ 25 % uncertainty should be taken into account by any analysis which either 
does not include final state interactions or cannot reproduce the numerical 
values of the phases. To reduce the uncertainties in the 1/N C approach it would 
be interesting to investigate the two- loop imaginary part and/or tcL-pambine 
our calculation with a dispersive calculation along the lines of refs.HcircJ. In 
order to reduce the scheme dependence in the result, appropriate subtractions 
would be necessary (see refs. Eilcj). Finally, it is reasonable to assume that 
the effect of the pseudoscalar mesons is the most important one. Nevertheless, 
incorporating the vector mesons and higher resonances would be desirable in 
order to improve the treatment of the intermediate region around the rho mass 
and to show explicitly that the large enhancement we find at low energy in the 
treatment of the pseudoscalars remains valid up to the scale ~ m c , where the 
matching with the short-distance part can be done more safely. 
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